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Abstract. In this paper, we study limit cycle bifurcations for a kind of non-smooth poly- 
nomial differential systems by perturbing a piecewise linear Hamiltonian system with a 
center at the origin and a homoclinic loop around the origin. By using the first Melnikov 
function of piecewise near-Hamiltonian systems, we give lower bounds of the maximal 
number of limit cycles in Hopf and homoclnic bifurcations, and derive an upper bound 
of the number of limit cycles that bifurcate from the periodic annulus between the center 
and the homoclinic loop up to the first order in e. In the case when the degree of perturb- 
ing terms is low, we obtain a precise result on the number of zeros of the first Melnikov 
function. 
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1 Introduction and main results 

There are many problems in mechanics, electrical engineering and the theory of au- 
tomatic control which are described by non-smooth systems, see for instance the works 
of A.F.Filippov [1], Andronov et al [2], Kunze [3] and the references therein. Recently, 
a good deal of work has been done to study bifurcations in non-smooth systems includ- 
ing Hopf, homoclinic and subharmonic bifurcations, etc.. In [4-6] Hopf bifurcation for 
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non-smooth systems was studied by developing new methods for computing Lyapunov 
constants. The Melnikov method for Hopf and homochnic bifurcations was extended to 
non-smooth systems in [7-9]. The method of averaging has also been extended to non- 
smooth systems in [10]. However, so far there are few papers in the literature studying 
homoclinic bifurcations inside the class of piecewise polynomial differential systems. In 
this work, we study this problem by using the first order Melnikov function of piecewise 
near-Hamiltonian systems. 

In [7], Liu and Han considered a general form of a piecewise near-Hamiltonian system 
on the plane 

I -T" — J-T _L cnt nr" -7/ A 1 



X = Hy + ep{x,y,6), 



where 



y = -H^ + eq{x,y,6), 



Tr(^ _ / H+{x,y), X > 0, 
^ I H-{x,y), x<0, 

p(x,y,d) = < ^ _) - „' 



q{x,y,6) 



q^ix,y,6), X > 0, 
q'ix,y,6), X < 0, 

H^, p^ and are C°°, e > is small, 6 E D G R"^ is a vector parameter with D 
compact. This system has two subsystems 



X = + ep-^{x,y,S), ^ , 



and 

X = Hy +ep~{x,y,5) 



y = -Ht + eq^{x,y,5), 

y = -H^ +eq [x,y,d), 

which are called the right subsystem and the left subsystem, respectively. Suppose that 
fll.ip l^-n has a family of periodic orbits around the origin and satisfies the following two 
assumptions. 

Assumption (I): There exist an interval J = and two points A{h) = (0,a(/i)) 

and Ai{h) = (0,ai(/i)) such that for h E J 

H+{A{h)) = H+{Ai{h)) = h, H-{A{h)) = H~{Ai{h)) = h, a{h) ^ ai{h). 

Assumption (II): The subsystem fll.lh )|c— n has an orbital arc starting from A{h) and 
ending at Ai{h) defined by H~^{x,y) = h,x > 0; the subsystem fll.lb )|^-n has an orbital 
arc starting from Ai{h) and ending at A{h) defined by H~{x,y) = H'^{Ai{h)),x < 0. 

Under assumptions (I) and (II), (l.l)|£=o has a family of non-smooth periodic orbits 
Lh = L'^ U Lf^, h E J. For definiteness, we assume the orbits Lh ioi h E J orientate 
anticlockwise. See Fig. 1. 

By Theorem 1.1 in [7], the first order Melnikov function of system (11. ip has the form 

^(^''^) = Tr77T W7t4 / q^dx-p+dy+ q'dx - p'dy , h E J. (1.2) 
Hy [A) lH+{Ai) Jl+ Jl- J 



Ai 




Fig. 1. The closed orbits of (l.l)|£=o- 



Also, we know from [7] that if M{h,6) has at most k zeros in h on the interval (a,/?) 
for all (5 G -D, then fll.ip has at most k limit cycles bifurcated from the open annulus 

Let 

M^{h, 5) = q^dx — p^dy, M''{h, 5) = q^dx — p^dy, 



M {h, 6) = ^ q dx — p dy, 

JH-(x,v)=h, x<0 



;i.3) 



'H-{x,y)=h, x<0 

where h is given in assumption (I). Then we know 

M-{h,6) = M-{H-{Ai{h)),6), he J. (1.4) 

As in the smooth case, a very important issue associated with (11. ip is to find the number 
of limit cycles and their distribution. In [7] , Liu and Han consider a piecewise polynomial 
system of the form 

(x y) = I ^^^y + ^' + ^' ^ - °' (1 5) 

1 {h~y + ep-{x,y,5),-h-x + eq-{x,y,5)), x < 0, 

where > 0, p^ and q^ are arbitrary polynomials of degree n. It was proved that 
the maximal number of limit cycles of the above system is n up to the first order in e. 
We know the unperturbed system of (II. 5p has a global center at the origin. Then, an 
interesting question is how many limit cycles system (11.11) can have if (11.11) |c—n is a linear 
piecewise system with a homoclinic loop and p^ and are arbitrary polynomials with 
degree n. 

In this paper, we take 

H+{x,y) = ^{{x-l)^-y^), x > 0, (1.6) 



H ix,y) = --{x^ + y^), x < 0, 



(1.7) 



and suppose 



p{x,y) 



q{x,y) 



p+{x,y) = (4jX'y\ X > 0, 

i+j=0 
n 

P'{x,y)= E arjx'y\x<0, 

i+j=0 
n 

q+{x,y) = J2 Kj^'y\ a: > 0, 

n 

(i'{x,y) = E Kj^'y\ x <o. 

i+j=0 



Under fll.6p and fll.7p . the system (11 .ip becomes 

X = -y + ep+{x,y), 
y = I - X + eq+{x,y), - ' 



X = -y + ep {x,y), 
y = x + eq-{x,y), 



X <0. 



;i.9) 



For system (11.91) |^-n. there exist a family of periodic orbits as follows 

Lh = {ix,y)\H-'ix,y) = ^, x > 0} U {(x, y) = ^, x < 0}, 

with h = h — l,0<h<l. For the sake of convenience, here we use h/2 instead of h. If 
h ^ 1, Lh approaches the origin which is an elementary center of parabolic-focus type (see 
[4] for the definition). And if /i — )■ 0, L/, — )■ Lq, where Lq is a compound homoclinic loop 
with a saddle 5*1(1,0). See Fig. 2. 




Fig. 2. Phase portrait of system (1.9)|£=o 



Note that Hy{0, y) = Hy {0,y) for —1 < y < 1. Then by (II. 2p . we have the first order 
Melnikov function of system (II. 9p satisfying 



^{~') ~ f q'^dx—p'^dy+ f q dx — p dy = M{h), 
2 JaAi JaTa 



:i.io) 



where < h < 1, and 



4 = (0, -VT^), A, = (0, VT^), 

AAi = {{x,y)\H-^{x,y) = l, x > 0}, 
A^ = {ix,y)\H-ix,y) = !^, x<0}. 

Let Z{n) denote the maximal number of zeros of the non-zero function M{h) on the 
open interval (0, 1) for all possible p and q satisfying (11 .SI) , which is the maximal number 
of limit cycles of (11.91) bifurcated from the periodic annulus Uo</i<i for all possible p 
and q satisfying (II. 8p when M{h) is not identically zero. Let Nnopfin) and Nhomociji) 
denote, respectively, the maximal number of limit cycles bifurcated in Hopf bifurcation 
near the origin and in homoclinic bifurcation near Lq for all possible p and q satisfying 
(II. 8p . Then our main results can be stated as follows. 

Theorem 1.1. For any n > 1 we have 

(I) NHopfin)>n+[^^]. 



(2) iVwoc(ri) >n+[^]. 



rn+il 



Theorem 1.2. For n = 1,2,3, A, we have Z(n) = n , ^ 2 j- 
Theorem 1.3. For any n > 5 we have n + [^] < Z{n) <2n+ [^]- 

2 Preliminary 

In this section we give an expression of the first order Melnikov function M{h) for 
< h < 1, together with the two expansions of M{h) near the origin and the homoclinic 
loop Lq, respectively. By (II. 3p and (II. 4p we have 

M+(-) = / q+dx-p+dy, M' {-) = / q' dx - p' dy = {^—) . (2.1) 



2 JaAi 2 2 



Then, applying Green's formula twice we obtain 



M^ib = - [ [ _ ___^ipt + q^)dxdy- [ p+{0,y)dy 

^ J Jint{AAi\JA^) JAAi 

= - f ,Pix,y)dy- / p^{0,y)dy, 

JAAilJAiA JAAi 



where 

pix,y) =p^{x,y) -p'^{0,y) + I q^{u,y)du, (2.2) 



satisfying 

p{0, y) = 0, p^=p+ + q+ 



Since AAi can be represented as x = 1 — y^h + y'^ = ilj{y, h), we have 

rVT-h 



4' 



^-{I,{h) + h{h)). 

n 

Note that p"'"(0, y) = Yl, ^ojV'' ■ follows that 



j=0 



p^{0,y)dy 

VT-h 



n + 

E7Ti<^-(-i)'"')(i-'')* 



k=0 



By (O, we have 



n n n 

i+j=0 j=0 i+j=0 



n-1 
i+j=0 

where 

In particular, for n = 2 we have 

Poo ~ '^10 + ^01' PlO ~ ^20 + 2^^"'^' ~ ~^ 2^02- 



The definition of Io{h) and fl2.5p yield 



P{i^{y,h),y)dy 

-VT-h 



t-VT-h 

i+j=0 J~VT^ 

2 E / (1 - ^^hTV'rv'dy. 

.: , n,._n "'0 



(2.7) 



For < A; < [^] and < r < n let 

Irk{h)= [^\h + y')h''dy. 
Jo 

Then for < A; < [^] and < i + 2k < n - 1 

•^0 r=0 

If r = 2Z, I > 0, we see 

hLk{h)= / {h + y^yy^Hy = VT^^i+k{h), (2.8) 







where ipi+k{h) is a polynomial of degree / + k. Hence l2i,k{h) is at /i = 0. If r 
2/ + 1, / > 0, then 



hi-,iAh)= [ {h + yj+h^'dy. (2.9) 
Using the formula 



we have from ( 12.91) that 

J ,(M = + (2^+1)^ J / > 1 A; > 0. 

21+iM J 2{l + k + l) 2{l + k + l) ^' - ' - 

It follows that 

hi+iAh) = (1 - h)''+-^^ik{h) + aikh'hkih), / > 1, > 0, (2.10) 

where 

(2/ + l)(2/-l)---3 



2^{l + k + l){l + k)---{k + 2)' 
and 

^^^^^^ 1 I (2/ + !)/^ , , (2/ + l)(2/-l)---5-/.'-i 



2(/ + A; + l) 2(/ + A; + l)(/ + fc) 2'(/ + A; + 1)(/ + A;) ■ ■ ■ (A; + 2) ^ 

which is a polynomial of degree I — 1 in h. For A; > 0, let 



* (h\ - ) Vik{h), / > 1, * _ / "ifc' ^ > 1' 

'^''^^''^ - ^ 0, / = 0, "''^ ~ 1 1, / = 0. 



By dllO]), 

/2m,fe(/^) = (1 - h)''+-'V*ikih) + aikh'hkih), l>0,k>0. (2.11) 



Further, by using the formula 
we have 

, , h(2k-l) , 

Thus 

hk{h) = VT^^k-i{h) + fSkh'hoih), k>l, (2.12) 

where 

^ , (2^-1)!! 

and 

(1-/^)^-1 (2fc - - Z^)^-^ , ,_, (2fc-l)(2fc-3)---3-/z^-^ 

^'-'^ ^ 2(A; + 1) 4(fc + l)fc + 1 2'=(fc + l)fc(fc-l)---2 ' 

which is a polynomial in h of degree k — 1. Let 
By (I2I2D, 

= VT^rk-iih) + Plh'^hoih), k>0. (2.13) 



Moreover, by (12.9 



r-y/l-h 2 2 

ho{h)= / (/i + i/^)^rft/ = -[Vl - /i + ln(l + Vl - /i) - -/I \nh] = — Mn/i + 7(/i), 
Jo 2 2 4 

(2.14) 

where 7(/i) is at /i = 0. We observe that 



ho{h) = h^o[\j^-f^), <^o(w) = ^ yr+^rfx. (2.15) 



The function '^^{u) is analytic on and odd in u. 
From (EH]), ( 12331) and ( l214ll . we can obtain 



= Vl - hifi+k-iih) + aikh^^'' Iio{h) (^2.16) 

= 7(/i) - la^,/i'+^+Hn/i, 

where i^i+k-iih) is a polynomial of degree I + k — 1, aik are non-zero constants, ^{h) is 
C"^ at /i = 0, < A; < [^], < 2Z + 1 < n. 



Substituting (EH) and fl216|) into ([2l]), and noting that k+[^] < [f ] and k+[^] < [^] 
since i + 2k < n — 1, we have 



n— 1 / i+l j+1 \ 

j+2fc=0 ^ r=0, r even r=l, r odd ^ 

n-1 



j+2fc=0 



(2.17) 



10 



= Vl — hu[^]{h) + z/jn--i|(/i)J 
= 0i(/i)ln/i + 02(/i), 

where i^j{h) denotes a polynomial of h of degree j, 02 (^) ^ C''^ at /i = 0, 0i(/i) is a 
polynomial in h of degree [^^] and 0i(O) = 0. 

Now we consider M~(|) in (12. II) . Set x = \/ —h cosO, y = V—h sinO for —l<h< 
0, ^ <e < and let 



2 — — 2 ' 



37r _ 
2 .,.12 



I.. = I cos'e sin^ede = {-ly^^ / cos^e sin-^'^rf^ = (-1)^+^1^^, 

where 



r 



0, j odd, 

a positive constant, j even. 



It follows from (11.31) that 



_ h 

M (— ) = / q dx — p dy 



2' 



n „ 32L 

i-h)'^ / ' (br^cos'e sin^+^e + a^.cos^+i^ sin^'^)d^ 



2+J=0 2 



i^ijkj+i + (2.18) 

i+i=0 



I 

i+i 



1=0 j=0 



where 

I 

3=0 

I 

= (-1)'+^ + (2.19) 

j=0 

I I 

j=0, j even i=l, i odd 



By dZH) and ([218]), 



2 

where for n = 2 



h " i 

M-(-) = -Vl - /^^^'(l - 0</i<l, (2.20) 



z=o 



7r__ 2___ 
eo = -2aoo, ei = -(a^o + ^01). ^2 = -g(2a^o - Ki + %2)- 

Hence by f lLTOj) . fl2J|) . flOj) . ([Ml), fl2Tri) and fl2:20|) we obtain the following lemma. 
Lemma 2.1. For system (II. 9p . the first order Melnikov function has the following form 

M{h) = yJT^fr,{VT^) + (7[^](1 - /i)/io(/i), 

where fn and g^n^^ are polynomials of degrees n and [^^] respectively. 

In the following we study the expansions of M{h) near h = 0,1. By (12.140 and Lemma 
2.1, the lemma below holds immediately. 

Lemma 2.2. For system ( \1.9\\ the first order Melnikov function M{h) has the following 
expansion near the homoclinic loop Lq 



M{h) = ( J2 b*ih') \nh + J2 bjh^^ < /i< 1, 

i=l j>0 

where for n = 2 

bo = — fli'o ~ ^01 ~ I ^20 ~ I ~ 2 flfj] — I a,^ + 2aQQ — |a^Q — ^Bq^ + |a2o — f&n + ^0-02^ 



2) + ^00 + ^02 '^00 



61 = (i + ln2) a+ + (I + ln2)5o+ + (21n2 - l)a+ + (ln2 

+ ^CLiQ + f &0I ~ 2a2o + ^11 ~ '^025 
"2 — 4 "20 "T 8 11 4 02 8 10 8 01 "T 4 "OO 4"00 2"02 4"!! "r 4"02) 



Noticing that in (I2.15P v'o(w) G on R and is odd in u, we can write for \u\ small 

oo 

where q is a constant, z > 0. Then, it follows from fl2.15p for 1 — /i > small 



ho{h) = h^o{^^) = = v^r^03(l - /^), (2.21) 

i=0 

where (psiu) G C"' at n = 0. Following (12.211) and Lemma 2.1, we have 

Lemma 2.3. For system 1^1. 9\} the first order Melnikov function M{h) has the following 
expansion near the origin 

n 

M{h) = yr^[^Ci(l -/i)5 +^c„+,(l -/i)I^]+^], < 1 - 1, 

i=0 j>l 

where for n = 2 

Co = 2(aoo - floo), ci = -f (aro + ^oi). 
C2 = -f (a^2 + ^01 + «w - 2a^o + - ^02), 
= -^(2^0 + «K) + + ^01). 

3 Proof of the main results 

Proof of Theorem 1.1. (1) For simplicity, we let and in (11.81) satisfy 

n 

p^{x,y) = ^afx\ q^{x,y) =0, (3.1) 

i=0 

that is, p'^i^x^y) are univariate polynomials of the variable x. By (12.41) and (12. Sp . we have 

n 

Ii{h) = 2a^Vl - h, p{x,y) = y^a+x\ 

From (II. 6p . the equation H^{x,y) = f is equivalent to 

2x-x'^ = u, (3.2) 

where uj = 1 — h — y"^. Consider equation (13. 2p in x near x = 0. There is a unique C°° 
solution 



i^ + 0(a;2)^^^y^ (3.3) 



X = a[uj) 2 



where ^ ^ ^ ^ ^ ^ ^ ^ 

Then, ( 13 ■3p imphes that the arc AAi can be represented by x = near the origin. 

That is, in this case the function ip{y, h) in (12. 3p can be taken as x = (t{u) with u = 
1- h - y'^, 1 - h> small Note that 

n 

i=i j>i 

where pi = ^af, p2 = ^a^ + ^af, • ■ ■ , pw = fi- + L{al, ■ ■ ■ , a+_ J, ■ ■ ■ , throughout this 
paper L(-) denotes a hnear combination. Then by (12.61) . we have for < 1 — /i ^ 1 



= 2J2pj [^\l-h-yydy 
i>i 



where i?j = 2 /q^(1 — u^Ydu is a positive constant for j > 1. Hence, by (12. 3p 



M\^) = -{h{h) + h{h)) 



-(Vl -h^PjBj{l - + 2a+Vl - /i) 



i>i 



(3.4) 



= Vl - /i^c*(l - /i)^ for < 1 - /i < 1, 



i>o 



where c* = <^ _~ r ■ >'i Under ([31]), we also have from fl2l9|l and fl2:20|l 
I J ^ -'-• 



M-(|) = -yr^^Q(l-/i)i 0</i<l, (3.5) 



where e; = (— l)'+^aj Ti+i.o, r«+i,o is a positive constant, < / < n. Thus, by (ll.lOp . (13.41) 
and ([33]) 

n 

M(/i) = v^r^[^Ci(l-/i)t + ^c„+j.(l-/i)[t]+J], 0<l-/i<l, (3.6) 

j=0 j>l 



with 

Co = Cq — Co = — 2a(j" + Fiocto j Ci = — ei = — r2o0.i , 



C2 = c[- 62 = — —aj + Tsoaa , C3 = -63 



C n 6^ 
2 



C„ 



-r4oa3 

l)"r„+i,oa; + L(a/ 
l)"r„+i,oa-, n odd, 



Br 
■ n 

2- 



2 



Cn+1 — C[n]_|_i — — 1 + 1 '^[fl + l ^'^1 ' ' " ' '^[f ]''' 



C2n-[f] = C„ 



Following the formulas above, we have 

5(co,Ci,C2, • • • ,C„,C„+i,--- ,C2n-[|]) 



where 









d{4 


1 ' ' ' ) '^[S] 


"[f]+i' 






■ ,an) 






-2 





■ ■ 











Tio 
















• • 














— r2o • • • 










Bi 
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• • 
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•• X* 

















(— i)"r„+i^o 








* 
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* 
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* — 




-^[f ]+2 


















v 
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Bn 

2" 












X* 



27 

0, n odd. 



n even, 



Clearly, the rank of this matrix is 2n — [^] + 1, which means that Cj, < i < 2n — [|] can 
be chosen as free parameters such that < |co| ^ |ci| ^ |c2| ^ • • ■ ^ |c„| ^ |c„+i| ^ 
|cn+2| ^ ■ • • ^ |c2„_[:|]| ^ 1, and QQ+i < 0, < 2 < 2n — [|] — 1 . Hence, by (13.61) . in 
this case M{h) has 2?2 — [|] positive simple zeros in (0, 1) near h = 1. Therefore, system 
(11. 9p can have 2n — [|] (=n + [^^], n > 1) limit cycles near the origin, which means 
NHopf{n)>n+m^]. 



(2) Here, we suppose and in (11.81) satisfy 

n n n 

p+{x,y) = ^aly\ q^{x,y) = ^bly\ p~{x,y) = ^a-x\ q~{x,y) = 0, (3.7) 



i=0 



1=0 



i=0 



where p^{x, y) and q^{x, y) are independent of x and p (x, y) is independent of y . Then, 
from ([MD to (1234]) . we have 



n-l 



p{x, y) = x^ bfiy' ^ = x^p^^y^ 

j=0 



i=l 



where 



Moreover, 



Poo = bt,Pm = • • ■ ,Ptn-i = ^K- 



k=0 



and 



3=0 



y'^)y^dy 



fc=0 



k=0 



'l-h 



k=0 



'l-h 



VhTWdy 



(3.8) 



fc=0 



fc=0 



where hk is given by Hence by ([O]), fl2T^ . (EH]) and the above 



E5lti(i-'")'-^E^(i-M'-E^oV:-.('") 



fe=0 



fc=0 



A:=0 



r "-l i 



(1 - h)-'^ Yl Phkf^lh'^hoih) 



k=0 



where 



n — 1 1 r n—l 1 



k=0 k=0 k=0 



\ 7l—l 

2 



k=0 k=0 



\ n—l 
2 



4 

k=0 



= ^^^^ + 5Z ^fc^'^^^ln/i, for < /i< 1, 



fc=0 A;=0 



Vq — Ojq -\- L[a2 , ■ ■ ■ , O^jnpP^lo, ■ • • ;P^2[ "-l ])' 
"^1 ~ ^ + -^^(,0-4 , ■ ■ ■ ) '^2[2.]'P00' ■ ■ ■ 'Po,2[^]''' 

^4: J_ T f + + + + \ 



V*2 = ^P04f^2+ L{P0Q,P02), 



* _ 1 



(3.9) 



^^()^2[Ii^]^Pfi] -^(^'do'^'d2' ■ ■ ■ '^'^2[^]-2''' (3.10) 

Moreover, from (13. 5p we have for < /i ^ 1 



(3.11) 



where 



Ci = I E 62.(1 - hy = -I E a2-,r2.+i,o(i - hy 

i=0 i=0 



(3.12) 



2C2 = E e2,+i(i - hy = E s+ir2j+2,o(i - hy 

j=0 j=0 

= VT^ E «2,+ir2,+2,o E c^i-i)'h'' (3.13) 

j=0 k=0 



Let 



where 



Denote 



v^T^ E ( E a2,+ir2,+2,oCj=(-i)'^)/i* 

A:=0 ^ j=ifc ^ 

Vl - /i = < /K 1, (3.14) 



j>0 



1 (2i-3)!! 
Ho = 1, /ii = --, /ij = , I > 2. (3.15) 



[^1 



Dk=Yl a2„.+ir2™+2.oC^^(-l)', 0<k< [\^]. (3.16) 



m=A; 

Then, by f l37[3|) . fl3:Til) and fl3:T6|) we have 



[^1 i [^1 



2C2 = 5Z (^/x.-,Z^,)/i' + E ( E /"fc+s^lV]-)^''^'""'' < ^ « 1- (3-17) 

i=0 j=0 k>l s=0 

From (l3ll]) . (I312|) . (I3A6|) and f lXTTj) . we get 

M~(-) 

-—^M= = J2^,h\ < « 1, (3.18) 
-2Vl-/i 

where 



2 J' 



/ n -l 

2 



C<7r n — 1 1 



I /Ui+s-Dn^_^, n odd, 



s=0 



71-2 

2 



K s=0 



r "-i i 



n-2 
2 



^'-^^Z r^+i,oa~, n even, 



(3.19) 



Hence, it follows from (13. 9 P and (I3.18P that 

J^i^ih' + v*h'+Hnh), < /K 1, 



M(|: 



(3.20) 



i>0 



where by (l3A0il and ([319]) 



3 -\- L[a'^ , ■ ■ ■ ,fl^p],Poo; 



'1 - A 



iPoVi* + ^(Poo), 



-^(-1) 2 ,P(|;„_i)+L(ao,ai,--- ,a„), nodd, 

^(-l)i-i + L(a+,p+,p+,--- ,p+„_2) + L(ao,a^,--- ,0"), n even, 



^ n -1 

2 



I E + ^(P00>P02, • • • ,P(tn-l), ^ odd, 

3=0 ^ 

|ij(-l)t+L(p+,p+,--- ,p^„_2) + L(a-,ar,aJ,--- ,a-_i), n even, 

0(|f*,f*,--- ,vr^A), i > 1, (obtained by dSH) and (13J0|) ). 

I 2 I 



s=0 



(3.21) 



If n is odd, by the formulas of f [nziij+j, j > 1 in the above and D^, < k < [^^] i 



m 



fl3.16p we have 

~ I ^ + + + 



n— 1 n — 1 

~2~ ~2~ 



J2 Yl fJ'j+sa2m+1^2'm,+2,oCm^ (-1) 2 + L{p'^Q, Pq^, ' ' ' , Po,n-l] 



n—1 n—1 
/ 



2 / 2 n—1 ^ n—1 \ 

m— n V . ri— 1 / 



2 

m=0 ^.,=Ii_^_m 



n-l 
2 



1 ^ /" ^ 

= 2 Z] ( /^j+IlzJ._fr2m+2,oCm(~l) ) '^2m+l + -^(PoO' Po25 ' ' ' ^^Cn-l); 
m=0 ^t=0 ^ 

(3.22) 

where j > 1. For 1 < j < it follows from ( KT2^ 

d(Vn-l,-. , Vn-1 , 9, ■ ■ ■ , Vn) 

^1 - ' ^ ^ " (3.23) 



n-l \ 
1 — \ 

2/^i+n^r2o 2 Z] /^i+ii^-tr4oCi(— 1) ••■ 2 Z] /^i+B^-^rn+i.oCn-i (— 1) 



1 

2/^2+21^^20 2 S /^2+ii^ -4^40^1 ( — 1) ••• 2 Z] /^2+S^-trn+l,oCn^( — 1) 
t=0 ^ t=0 ^ 2 



n-l 
2 



1 



n-l 
2 



|Ai?ir2o I X] /^n-t-'^4oC'i(— 1)* ••• I ^ /i„_ir„+ioCi_^(— 1) 



\ t=0 t=0 2 / 

Let Q denote the ith column of Ai, a x q (or q/o) indicate multiplying (or dividing) 
each element of the ith column by a, and q + 6 x cj indicate that each element of the jth 
column times b adds to the according element of the ith column, a and b are constants. 
We make elementary transformations to Ai with the following steps 

51. Q/(ir2.,o),2 = l,2,---,^. 

52. c, + (-l)xci,z = 2,3,---,^. 

53. Ci + i-Cl,) X C2, z = 3, 4, ■ ■ ■ , (-1) x c^. 



Sn + l . C n+1 -|- ( — ^^) X Cn-1, Cn_J,/( — 1) 2^ aud Cn±l/( — 1) 2 , 



Then Ai becomes 



Ai 



n-l 



/i2 /il \ 

yUs /^2 



(3.24) 



By f l3.15p . we can prove Rank(74i) = See Appendix for the proof. Therefore, by 

using f l3.2ip and f l3.24p we obtain 



(9(ag , a2 , ■ ■ ■ , C^ri- 15^00)^02 5 ■ ■ ■ 5Po,n-l5 '^1 5 '^3 ' ' ' ' 5 '^n ) '^0 ' '^2 5 ' ' ' 5 '^n-l) 



/ 


-1 

1 


* 


* 


* 


* 


* 


* 


* 


* 


* 


* 


* 


\ 




n 
U 


1 

3 


* 


* 


* 


* 


* 


* • • 


* 


* 


* 


* 










n-1 



























•• 


(-1)^ 

n 


* 


* 


* 


* 


* 


* 


* 


* 


* 









•• 





1 

4 


• 








•• 


• 





• 


• 









•• 





* 


4 








•• 


• 





• 


• 





































•• 





* 


* 


4 





•• 


• 





• 


• 









•• 





* 


* 


* 











• 


• 
































v 





•• 





* 


* 


* 











• 


• 


/ 



It is easy to see that the rank of this matrix is n + + 2. Thus, we can choose 
vo,Vq,vi,vI, ■ ■ ■ . Vn-1 . f *-i , Vn-1 I ^ ■ ■ ■ ,Vn a.s free parameters such that 

2 2 2 

< Wo < ^^n < ^'i < < ■ • • < < ^"-1 

2 — 2~ 

, , , , . . n — 1 , 

< Vn^+l < ~Vn-i_^_2 < ■ • • < (-1) 2 t;„ < 1, 

or 

" 2 —2" 

^ , . , ^ n-|-l ^ , 

< -■Wb-J.+I < < ■ ■ • < (-1) 2 < 1- 

Note that in (jMI]) 

. = 0{\vl,v\, ■ ■ ■ I), i > 1. 

2 2 

Then, by fl3:20|) M(/i) has n + ^ positive simple zeros near h = 0. Hence, system (11. ip 
can have n + (= n + [^^] with n odd) hmit cycles near the homoclinic loop Lq. 

If n is even, by using the same method as above, we get the rank of the following matrix 

d{aQ , a^, ■ ■ ■ , a+, p^, Pq2, ■ ■ ■ , '^1 5 '^3 ' ' ' ' ' ^n-i^ '^0 ' '^2 > ' ' ' j^n) 
is n + 1 + ^, which means this matrix has a full rank in row. Thus, we can choose 

< < t^n < ^1 < ^'t < ■ ■ ■ < ^li^ < ^* -2 

2 

< t;^ < -vn+i < • ■ ■ < (-1) < 1, 

or 

< -vo<^ -vX < -1^1 < -t;? < ■ ■ • < -t^n^ < -t;* -2 
" 2 



Also by (K20\i and vl_2^. = 0{\v^,vl, ■ ■ ■ ,<_2|), i > 1, the function M{h) has n + ^ 

2 2 

(= n + [-^^J with n even) positive simple zeros satisfying < /i ^ 1. Therefore, system 
f ll.Qp can also have n + [^^] limit cycles near the homoclinic loop Lq in this case. That 
is, Nhomociji) >n+ [^^]. The proof ends. 

Proof of Theorem 1.2. By Theorem 1.1, it suffices to prove that Z{n) < n + [^^] 
for = 1, 2, 3 and 4. 

Let n = 1. Then by Lemma 2.1 and fl2.15p 



M{h) = AVl - h + 5(1 -h) + Chipo 
where A,B,C are constants. We have 



h 



M{h) = Vl - h[A + BX + C^i(A)], 



where A = y/1 — h G (0, 1), and v^i(A) 



A 



. Note that 



A 



1 - A 



2^2 



It is direct that 



and 



Note that 



IM^)]' = -{1 + x-')vo 



A 



A 



+ 



A(l-A2)' 
A , 



1-A2 



u 



\/l + u"^ if u 



A 



We have v^o(^) < mVT+I? since {(foiu) 



u\/\ + v?y < for M > 0. Therefore [y9i(A)]" < 0, which implies that M{h) has at most 
two simple zeros in the interval (0, 1). That is, Z{1) < 2. 

For the case of n = 2, it follows from Lemma 2.1 and (12.151) that 



M{h) = Vl^[Ai + A2\ + AsX^ + A^ifiiX)], 



where A = Vl — h G (0, 1), Ai, A2, A3 and A^ are constants. By simple computation, we 
get 

A(4A2-3)- 



6 



smce 



A 



+ 



A 

A(4A2 - 3) 



+ 



3(1- A2)2 
1 

(1-A2)2 



< 0, 



> 0. 



~VT^^ 3(1-A2)2. 
Hence, M{h) has at most three simple zeros in the interval (0,1). This means Z(2) < 3. 
For = 3, we have 

M{h) = VT^IAq + AiX + A2A' + AsX' + {Bo + BiX^)^i{X)] = VT^Mi{X). 



Hence, 



We can find 



and 



)(A) = 12B^if'|{X) + 8B,Xif'l'{X) + {Bo + BiX'')if^^\X) 
= B,{l2ip'l{X) + SXip'I'iX) + AVS'^(A)) + Boip^^\X). 

X \ A(5A^ - 7X^ + 3) 



VT^^ 3(1 - A2)3 



wax) + 8A<(A) + AVf ^(A) = 



(1 - x^y 



Thus, 



Since 



((l-A^)Vf)(A))' = 



(1- A2)3 

24(l-A2)2(A2 + 5) 



A6 



A 



(1 - A2)3 

A(17A2 - 15) 



and 



X \ A(17A2 - 15) 



+ 



y/T^J 3(1-A2)2(A2 + 5)J 

> 0, (3.25) 



.^/r^p; 3(1-A2)2(A2 + 5)J (1-A2)3(5 + A2)2 
we have for A e (0, 1) 

M^(A) 7^ if So 7^ 0. 

Therefore, M{h) has at most 5 zeros in the interval (0,1) for this case. It follows that 
Z{3) < 5. 

For n = 4, by using the same method above, we have 



and 



where 



and 



M{h) = ^/^^ ^Ao + A + ^2 A2 + A^X^ + AA^ + {Bo + BiX^)ipi{X)\ 
= VT^M3(A), 

(A) = 20Bi<p'l'{X) + 10BiX<p^^\X) + {Bo + BiX'')<pf\X) 
= i3i(20(^i"(A) + 10Xipf{X) + AVP(A)) + Boipf{X), 

A(30A6 - 59A^ + 50A2 - 15) 



(3.26) 



2099;"(A) + ioA<^f)(A) + aV^Ha) 



Thus, 



Mf (A) 



A 

(1-A2)^ 



- 48Si + Bo 



15(1 - A2)4 

48A 
~ ~(1-A2)4- 

^^^^f^A) 



(3.27) 



Noting that 



where 



we have 



since 



'1-A 



2\4 



A 



2\2 



960(7 -A2)(l-A2) 



M4(A), 



M4(A) = ipo 



A 



r (i-Ay 

. A 

M4(0) = 0, M^(A) 



+ 



A9 

A(48A^ - 155A^ + 105) 
15(A2-7)(1 - A2)2 



< 



(3.28) 



IGA^ri - A^ 



> 0. 



5(1-A2)4(7-A2)2 

Therefore, by (13.261) . (I3.27P and (I3.28p . M{h) has at most 6 simple zeros in the interval 
(0,1), which yields Z{A) < 6. The proof ends. 

Proof of Theorem 1.3. Following Theorem 1.1, we only need to prove Z{n) < 
2n + [^^], n > 5. We give the proof using the same method as in [11]. Let A = 
Vl — h, h G (0, 1). Then by (I2.15P and Lemma 2.1, it follows that 



M{h) = VT^UVT^) + (7[^](l - h)ho{h) 

= A/„(A) + (1 - A2)^[^,(A2)^o(7i^) 

= A/„(A) + n(A)^o(^) 
= M*(A), 

where u{\) = (1 — A^)mo, uq = (7[iLzi](A^). Hence, 



d .M*(A). 



A/„(A) 



u 



with 



^9(A) = {Xfn)'u - XfnU' + Ul 



(3.29) 



By fl3.29p and degree(M)=2[^^], we find that the coefficient of the term of ^{X) of degree 
n + 2[-^^] vanishes if n is odd. Thus, degree ('(9)< 'n, + 2[-2Yi] — 1 if n is odd, and degree ('i9)< 
n + 2[^^] if n is even. That is, degree('i9)< 2n. Let / = (0,1). Using notation jj{A G 
J|/(A) = 0} to indicate the number of zeros of the function / in the interval / taking into 
account their multiplicities, we have 



tl{AG/|M(A)=0}<[ 



n — 1, 



tj{A G I\^{X) = 0} < 2n. 



Therefore, 

tl{A G /|A^*(A) = 0} < G I\u = 0} + G /|F = 0, m ^ 0} + 1 

< tl{A G I\u = 0} + tl{A G /|i9(A) = 0} + 1 

,n + 1, 
<2n+[^]. 

The proof is completed. 



Appendix 

The proof of Rank(y4i) = (for n odd). For simphcity, we take n = 7 for 
example, since the method used here can be generalized for any n > 1. Thus, by fl3.15p 
and ^Ml, 

( -5!!/24 -3!!/23 -1/22 _i/2 \ 
r _ -7!!/25 -5!!/24 -3!!/23 -1/22 
^~ -9!!/26 -7!!/25 -5!!/24 -3!!/23 " 

\ -ll!!/27 -9!!/26 -7\\/2^ -5!!/24 / 

Let rj(cj) denote the zth row(column) of Ai for i = 1,2,3,4. Now we making elementary 
transformations to Ai. Multiplying rj by —2* for i = 1, 2, 3, 4, and cj by 2^"-' for j = 1, 2, 3, 
gives 



~ Ai 



( 5!! 3!! 1 1 \ 
7!! 5!! 3!! 1 
9!! 7!! 5!! 3!! 
\ 11!! 9!! 7!! 5!! / 

Further, for Ax, by adding rj times — (2z — 1) to rj+i for z = 3,2, 1, and dividing Cj by 
2(4 — j) for j = 1, 2, 3, we have 

/ * * * 1 \ 
5!! 3!! 1 
7!! 5!! 3!! 
\ 9!! 7!! 5!! / 

Using the same method as above, we also have 

/* * *l\ /* **l\ 

5!! 3!! 1 * * 1 

7!! 5!! 3!! ^ 5!! 3!! 

\ 9!! 7!! 5!! / \ 7!! 5!! / 



( 



* 1 \ 

1 



1 
\ 5!! / 



Therefore, Rank(y4i) = 4. The proof ends. 



References 



A.F. Filippov. Differential equations with discontinuous righthand sides. Kluwer Academic. 
Netherlands, 1988. 

A. A. Andronov, S.E. Khaikin, A. A. Vitt. Theory of Oscillators. Pergamon Press, Oxford, 
1965. 

M. Kunze. Non-smooth Dynamical Systems. Spinger-Verlag, Berlin, 2000. 

B. Coll, A. Gasull, R. Prohens. Degenerate Hopf bifurcations in discontinuous planar sys- 
tems. J. Math. Anal. Appl., 253 (2001) 671-690. 

A. Gasull, J. Torregrosa. Center-focus problem for discontinuous planar differential equa- 
tions. Int. J. Bifur. Chaos, 13 (2003) 1755-1765. 

X. Chen, Z. Du. Limit cycles bifurcate from centers of discontinuous quadratic systems. 
Comput. Math. Appl., 59 (2010) 3836-3848. 

X. Liu, M. Han. Bifurcation of limit cycles by perturbing piecewise Hamiltonian systems. 
Internat. J. Bifur. Chaos Appl. Sci. Engrg., 5 (2010) 1-12. 

Z. Du, W. Zhang. Mclnikov method for homoclinic bifurcation in nonlinear impact oscilla- 
tors. Comput. Math. Appl., 50 (2005) 445-458. 

F. Battelli, M. Feckan. Bifurcation and chaos near sliding homoclinics. J. Differential Equa- 
tions, 248 (2010) 2227-2262. 

J. Llibre, A. Makhlonf. Bifurcation of limit cycles from a two-dimensional center inside i?". 
Nonlinear Analysis, 72 (2010) 1387-1392. 

J. Llibre, H. Wu, J. Yu. Linear estimate for the number of limit cycles of a perturbed cubic 
polynomial differential system. Nonlinear Analysis, 70 (2009) 419-432. 

M. Han. On Hopf cyclicity of planar systems. J. Math. Anal. Appl., 245 (2000) 404-422. 

M. Han, G. Chen, C. Sun. On the number of limit cycles in near-Hamiltonian polynomial 
systems. Internat. J. Bifur. Chaos Appl. Sci. Engrg., 6 (2007) 2033-2047. 

Y. Wu, Y. Cao, M. Han. Bifurcations of limit cycles in a zs-equivariant quartic planar 
vector field. Chaos, Solitons & Fractals, 38 (2008) 1177-1186. 

J. Yang, M. Han. Limit cycle bifurcations of some Lienard systems with a cuspidal loop 
and a homoclinic loop. Chaos, Solitons &; Fractals, 44 (2011) 269-289. 

M. Han, W. Zhang. On Hopf bifurcation in nonsmooth planar systems. J. Differential 
Equations, 248 (2010) 2399-2416. 

M. Han, J. Chen. On the number of limit cycles in double homoclinic bifurcations. Sci. 
China, Ser A (2000)914-928. 

M. Han. Cyclicity of planar homoclinic loops and quadratic integrable systems. Sci. China, 
Ser. A (1997)1247-1258. 



